Astronomy & Astrophysics manuscript no. TanGoI'astro-ph 


©ESQ 2010 


January 27, 2010 





o 

(N 
G 

(N 



o 

u 

6 



> 
OS 

o 
o 
in 



TANGO I: 

ISM in nearby radio galaxies. Molecular gas 

B. Ocana Flaquer 1 , S. Leon 21 , F. Combes 3 , and J. Lim 4 ' 5 

1 Instituto de Radio Astronomia Milimetrica (IRAM), Av. Divina Pastora, 7, Nucleo Central, 18012 Granada, Spain 
e-mail: ocana@iram.es 

2 Joint Alma Observatory, Av. El Golf 40, Piso 18, Las Condes, Santiago, Chile 

3 Observatoire de Paris, LERMA, 61 Av. de l'Observatoire, F-75014. Paris, France 
e-mail: francoise . combe s@obspm. fr 

4 Institute of Astronomy and Astrophysics, Academia Sinica,128, Sec. 2, Yen Geo Yuan Road, Nankang, Taipei, Taiwan, R.O.C. 
e-mail: lim@asiaa. sinica . edu. tw 

5 Department of Physics, University of Hong Kong, Pokfulam Road, Hong Kong 
e-mail: jjlim@hku.hk 

Received MONTH day, YEAR; accepted MONTH day, YEAR 

ABSTRACT 

Context. Powerful radio- AGN are hosted by massive elliptical galaxies which are usually very poor in molecular gas. Nevertheless 
gas is needed in the very center to feed the nuclear activity. 

Aims. We aim to study the molecular gas properties (mass, kinematics, distribution, origin) in such objects, and to compare them with 
results of other known samples. 

Methods. We have performed at the IRAM-30m telescope a survey of the CO(1-0) and CO(2-l) emission in the most powerful radio 
galaxies of the Local Universe, selected only on the basis of their radio continuum fluxes. 

Results. The main result of our survey is the very low content in molecular gas of such galaxies compared to spiral or FIR-selected 
galaxies. The median value of the molecular gas mass, including detections and upper limits, is 2.2 x 10 s M . If separated into FR-I 
and FR-II types, a difference in H 2 masses between them is found. The median value of FR-I galaxies is about 1.9 x 10 s and higher 
for FR-II galaxies, about 4.5 x 10 8 M but this is very probably entirely due to a Malmquist bias. Our results contrast with previous 
surveys, mainly selected through the FIR emission, implying larger observed masses of molecular gas. Moreover, the shape of CO 
spectra suggest the presence of a central molecular gas disk in 30% of these radio galaxies, a lower rate than in other active galaxy 
samples. 

Conclusions. We found a low level of molecular gas in our sample of radio-selected AGNs, indicating that galaxies do not need much 
molecular gas to host an AGN. The presence of a molecular gas disk in some galaxies and the large range of molecular gas mass may 
indicate different origins of the gas that we can not discard at present: minor/major merger, stellar mass loss or accretion. 

Key words. Galaxies: evolution - Galaxies: luminosity function, mass function - Radio continuum: galaxies 



1. Introduction origin of the molecular gas. 

. It is widely believed that the AGN are powered b y accretion 

O . Radio galaxies have strong radio sources in the range of 10 41 f ISM onto a supe r massive black holes (SMBHs) (Antonucci 

1 1/yt6 —.„/„ »nn;,7olonf t« IflM _ 1 »-J v v 3 



to 10 4b erg/s, equivalent to 10^ 4 - 10^ y Watts, normally hosted 1993) According to lde Ruiter et al.1 (120021) the presence of ISM 

by giant elliptical galaxies with visual luminosities of about fa the circumnuclea r regions of AGN is indeed inevitable and 

- _ , 2.1xl0 10 /T 2 L dKellermann & VerschuH 0988^, calculated the large . S cale dust/gas systems should be related to nuclear 

^ with a Hubble constant of H=100 km.s J Mpc ). Although activity 



H i in general relatively deficient in cold gas, elliptical galaxies More recent studies bv | Ha rdcastle etlil (120071) suggest that it is 

y* , still contain some dense and cold Interstellar Medium (ISM) posfjible that ^ the apparent i y different accretion modes may 

detected in CO(1-0) emission (e.g. | Wiklind & Rydbeck| | 1986D . be the result of a different sou rce for the accreting gas. This 

Dust is also detected in the far-infrared (FIR) radiation which hypothes i s came in lAllenet a l. (2006) who showed that some 

is assumed to be thermal and to originate from dust heated by low luminosity radio galaxies in the center of clusters could be 

young massivestars or by the active galactic nuclei (AGN) powered by Bondi accretion, a spherical accretion of the hot, 

(Knaooetal. 1989fc | Wiklind & Henkel| [1995|). According to x . ray emitting me dium. The supporting arguments about the 

l Kennicutl l U222D> early-type galaxies show no independent nature of the accl - e ti on mode i n low power radio sources were 

evidence of high star formation rates (SFRs), suggesting that discussed as well by Best et al.1 d2006h . 

the older stars or the AGN are responsible for much of the FIR For i ow .i uminosity AGN, the stellar component might be 

emission. sufficient to fuel the active nucleus, but for powerful radio 

I Wrklind&Henkell dl995|) show that in elliptical galaxies the galaxiefj large scale dynamical processes are required to transfer 

gas is unrelated to the stellar populations, and favor an external lar momenm m f the inter stellar gas across the disk radius 

( Shlosman eTaTIl 1 9901: ICombes 2002). 
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Fueling is also dependent on the morphological type of 
the galaxy, since the radial gas flows are driven by non- 
axisymmetric instabilities, and the stability of the disk is essen- 
tially depend ent on the bulge-to-disk ratio ( Mih os & Hernquisl 
1996, 119941) . In early type galaxies, for example, once the 
accreted gas has settled in a symmetric disk, the large-mass 
black ho les (BH) are st arving, unless perturbed by an external 
tide (e.g. lCombesll2002l) . 

In order to better understand the role of the Interstellar Medium 
(ISM) in the radio galaxies we have built two samples, at low 
and medium redshift, of radio galaxies selected only on the 
basis of their radio continuum emission. 

This is the first paper of the Thorough ANalysis of radio- 
Galaxies Observation (TANGO) project of a series of papers 
that will cover a wide range in wavelength (optical, radio, IR, 
X-Ray). 

An observational criterion of the TANGO samples at low and 
medium redshift is that most of the objects are visible from 
the main (radio) telescopes of the the Northern and Southern 
hemisphere, such as ALMA, IRAM telescopes, VLA, VLT, 
GTC among others. 

The primary issue addressed in the present paper is the study of 
the molecular gas through the CO emission in the low redshift 
TANGO sample of radio galaxies (z<0.1). We want to test, 
in particular, whether samples selected on the FIR emission, 
e.g. lEvans et all d2005l) . and radio continuum-selected samples 
behave in similar ways. 

This paper is structured as follows: we describe the sample 
of radio galaxies, their observations and their data reduction 
process in § [2] The molecular gas properties are described in § 
[3] while the dust properties are described in §[4] and the gas and 
dust relations are discussed in § [5] Section § [6] is dedicated to 
the comparison of the result of this work with similar samples. 
Finally the discussion and the conclusions are drawn in §[7] and 
§[8] respectively. Properties of individual galaxies are presented 
in the Appendix [A] the beam/source coupling is explained in 
details in the Appendix iBl and the Appendix ICl is dedicated to 
the detail analysis of one of the galaxies in TANGO, 3CR 3 1 . 



2. Observation and data reduction 

2.1. Sample 

As previously mentioned in the introduction, the galaxies se- 
lected in the low and medium redshift TANGO sample have one 
main criterium of selection: a high radio continuum power. Thus 
in the low redshift sample studied in this paper, all the radio 
galaxies have a radio power at 1 .4 GHz larger than 10 22 5 W.Hz 1 
with a median of 10 24 4 W.Hz 1 . This makes the main difference 
with the previous samples, chosen e.g. by their FIR flux (e.g., 
Evan s et al .2001. 

This work studies a total of 52 nearby radio galaxi es, most of 
them from the Thir d Cambridge Catalog (3CR - Edg e et al. 
19591: iBennettl fl962h . New General Catalog (NGC - iDrevei 
1888 ), and the Second Bologna Catalog of Rad io Sources (B2 
- ICollaetal.ll 197(1 1 19721 119731: iFanti et al.lfl974h . There is one 
galaxy from the Ohio State University Radio survey Catalog 
(OQ - IScheer&Krauslll967h and one g alaxy from the Uppsala 
General Catalogue of Galaxies (UGC - lNilsonlll973l) . The de- 
tails of the catalogs for the sample are summarized in Table Q] 





it ill' I -'lI'lXll'L' 


G£i ill' fliitiw'liiin J 
/( Ul UCICC11UII 


3CR 


24 


58% 


NGC 


13 


62% 


B2 


13 


31% 


UGC 


1 


0.0% 


OQ 


1 


100.0% 


Total 


52 


52% 



Table 1. The catalogs are listed in order of importance, some 
galaxies are in more than one catalog with different names. In 
this article the names are chosen first in favor of the 3CR cata- 
log and last for the OQ catalog. 1 This % is within the catalog, 
meaning the percentage of galaxies detected in each catalog. 



2.2. Observations 

This work is based on a survey of the 12 CO(1-0) and 12 CO(2-l) 
transitions, simultaneously observed with the IRAM 30 m tele- 
scope, at Pico Veleta - Spain, reaching an rms temperature of 
about 0.6 mK with a velocity resolution between 20 km/s and 52 
km/s. 

The transitions were observed at 115.27 GHz and 230.54 GHz 
(where the beam of the IRAM-30m telescope is about 22 and 1 1 
arcsec), redshifted at the velocity of the galaxies, for the 12 CO(l- 
0) and 12 CO(2-l) respectively with the A and B receivers. The 
filter backends were configured into two units of 512 x 1 MHz 
channels and two units of 512 x 4 MHz, one for each receiver. 
During the observations, the pointing and focus were monitored 
by observing planets and standard continuum sources every hour 
with an accurancy of about 3 arcseconds. 
From the total 52 galaxies in the sample, all of them were ob- 
served in the CO(1-0) line and 43 galaxies have simultaneously 
been observed in CO(2-l) line, when the weather permitted. 
About 58% of galaxies (N=30) in the sample have been de- 
tected (clearly or tentatively) in one or both lines and 38% of the 
galaxies (N=20) have been clearly detected, some of them only 
in CO(1-0), some only in CO(2-l) and some in both CO transi- 
tions. A total of 7 galaxies have been detected in both lines, these 
galaxies are 3CR 31, 3CR 264, NGC 326, NGC 4278, NGC 
7052, B2 01 16+31 and B2 0648+27. These 7 galaxies have been 
used to study the CO(2-l)-to-CO(l-0) line ratios that will be dis- 
cussed in section § I3.2l in more detail. 

There were 9 galaxies for which we only have data in the 3 mm 
wavelength range (3CR 236, 3CR 305, 3CR 321, 3CR 327, 3CR 
403, 3CR 433, NGC 6251, B2 0836+29 and OQ 208) where 6 
of them have shown a clear detection (3CR 305, 3CR 321, 3CR 
327, 3CR 403, B2 0836+29 and OQ208). 
Table|2]is a journal of the observations which shows the dates on 
which each detected galaxy was observed. 
This project started in 1999 and the last observations are from 
the year 2007, some of the galaxies have been observed sev- 
eral times throughout these 8 year period to confirm some pos- 
sible detection s. Part of this work has already been published by 
Lim et al ] (l2000h with a detailed study of the galaxies 3C3 1 and 
3C264. 

We observed as well the CO(1-0) emission in the center of 3C31 
using the IRAM Plateau de Bure Interferometer (PdBI). We used 
several tracks in the B configuration. The velocity resolution was 
set to 20 km s'and the CLEANed maps have been restored with 
a synthesized beam of 2. "2 x l."2 (PA=47°). See AppendixO 
for a clear analysis. 
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Galaxy 


Observation Dates 


Galaxy 


Observation Dates 


3C31 


1999 Dec 


NGC315 


2002 Jun 




2000 Jan 




2006 Oct 


3C6666 


2000 Jan 




2007 Apr 




nn o a it 

2003 May 


/V&C326 


1AA/" /" A . 

2006 Oct 


3C83.1 


2000 Jan 


/VG/C541 


2002 Jun 


3C88 


1999 Dec 


/V&C708 


2006 Oct 




2000 Jan 


/V&C3801 


2002 Jun 


3C129 


2000 Jan 


/VGC4278 


2007 Jan 


3C264 


2000 Jan 




2007 Apr 


3C272.1 


2000 Jan 


/VG/C5127 


2002 Jun 


3C274 


2000 Jan 


/V&C7052 


2002 Jun 


3C305 


2004 Feb 




2006 Jul 




2004 Mar 


62 01 16 + 31 


1AA1 r-\ 

2003 Dec 


3C321 


2004 Mar 




2006 Jun 


3C327 


1Art A ~\ It 

2004 Mar 




2006 Oct 


3C353 


2000 Jan 


62 0648 + 27 


2004 Feb 


3C386 


Yyyy Dec 




ZUU4 Mar 




900^ Dpp 




Z.UUU Utl 




2004 Feb 


52 0836 + 295 


2004 Mar 


3C442 


2000 Jan 


B2 0924 + 30 


2007 Apr 


3C449 


1999 Dec 


62 1347 +28 


2007 Jan 




2000 Jan 




2007 Apr 






O£208 


2004 Mar 



Table 2. 

Observations journal of the 30 detected galaxies in the sample. 



2.3. Data reduction 

For all the data reduction, the Continuum and Line Analysis 
Single-dish Software (CLASS) was used. All spectra were in- 
dividually checked, the bad ones were removed, and the remain- 
ing spectra of the same frequency were added. The baseline was 
removed after selecting a window. This window was taken to 
be the velocity width for the detected galaxies, and a line-width 
value of ~300 km/s was applied to compute the upper limits. 
Each spectra was smoothed to a velocity resolution between 20 
and 50 km/s. The integrated intensity, the velocity width, po- 
sition and peak temperature of the galaxies were determined by 
fitting a Gaussian. The result spectra of each detected galaxy can 
be found in the online material. 

We use the baseline of the spectra to calculate the continuum 
flux at 3 and 1 mm. The continuum level and the detections were 
given in antenna temperature, T*, and then converted to main 
beam tem perature, Tub- The con version factor used was the one 
given bv lRohlfs & Wilson! d2004t) . where we divided the forward 
efficiency (which is a model function to the measure antenna 
temperature and the copper load temperature) by the beam effi- 
ciency (which is the percentage of all power received that enters 
the main beam), F e ff/B e ff. Applying this for each frequency in 
the case of the IRAM-30m telescope, the relationship is 1 .27 at 
3mm and 1.75 at 1mm. 

To convert the data from temperature (K) to flux (Jy), the point 
source sensitivity measurement (S/FJ) is 6.3 for 3mm and 8 for 
1mm wavelength. The results for these conversi ons are listed i n 
Table [3] For more details, on the conversion see lKramerldl997h . 
From the observations, in summary, 90% out of the 52 galaxies 
were detected in the continuum at 3mm and 65% out of the 43 
galaxies were detected in the continuum at 1mm. 



2.4. IRAS Data 

We use the IRAS data to get the dust properties both for the 
warm component (A ~ 60/lotj) associated with young star form- 



Galaxy 


*J 3mm 


°^ 3fflffl 


c 

" 1mm 


rr 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


3C31 


79.07 


1.13 


74.39 


2.78 


3C40 


96.91 


0.9 


20.11 


1.8 


3C666 


110.12 


0.43 


73.08 


1.39 


3C83.1 


33.45 


0.69 






3C88 


121.84 


0.43 


87.05 


1.3 


3C98 


18.09 


0.52 


6.18 


1.54 


3C129 


28.53 


0.65 


0.91 


1.27 


3C236 


225.06 


0.57 


X 


X 


3C264 


174.64 


0.5 


87.35 


1.13 


3C270 


326.34 


0.82 


204.71 


1.57 


3C272.1 


138.29 


1.2 


73.86 


6.79 


3C274 


n a <i 1 o 

2347.13 


1.83 


1312.31 


2.78 


3C296 


99.23 


0.6 


54.55 


1.22 


3C305 


21.23 


4.28 


X 


X 


3C321 


67.33 


0.52 


X 


X 


3C327 


105.81 


0.69 


X 


X 


3C353 






3.46 


1.35 


3C386 


50.21 


0.88 


105.01 


2.26 


3C402 


332.57 


39.27 






3C403 


20.35 


0.48 


X 


X 


3C433 


91.29 


0.76 


X 


X 


3C442 


9.62 


0.43 


75.51 


9.99 


3C449 


25.82 


0.5 


15.83 


1.65 


3C465 


86.06 


0.56 


27.32 


1.31 


A If 1 1 C 

NGC3I5 


353.56 


0.67 


62.12 


1.39 


A I S~°* ST) 1 f 

/V&C326 


63.47 


0.63 






A T f~1 f~~< E A 1 

NGC54I 


4.73 


2.02 


27.23 


2.60 


/V&C708 


8.11 


0.98 


148.94 


3.83 


A TS~~1 S~"~\ A O A 

/V&C2484 


110.63 


0.69 


11.92 


1.65 


NGC2H92 


190.05 


0.87 


272.39 


2.23 


/V&C3801 


1 181.72 


210 






/V&C4278 


215.9 


1.08 


194.34 


1.28 


A TS~~1 S~~< E i n 

NGC5121 


25.62 


1.13 






/V&C5141 


63.50 


1.01 






A T f~i f~~< E A f\C\ 

/V&C5490 


27.72 


1.44 


88.61 


5.46 


NGC625 1 


367.88 


0.71 


X 


X 


/V&C7052 


85.51 


12.13 


18.77 


1.49 


B2 0034 + 25 


72.97 


0.68 


190.81 


1.82 


B2 0116 + 31 


164.53 


1.22 






B2 0648 + 27 


152.99 


0.37 


172.99 


2.09 


B2 0836 + 296 






X 


X 


B2 0915 + 326 


11.39 


1.29 


24.41 


2.17 


62 0924 + 30 


9.44 


1.13 






j > -\ 1 1 A1 1 O 

62 1 101 + 38 


359.89 


0.97 


0.165 


0.002 


62 1347 + 28 


8.96 


0.51 






62 1357 + 28 










do 144.7 _i_ 97 




U.J J 


S? 41 


n 81 

U.o 1 


62 1512 + 30 










62 1525 + 29 


14.48 


1.07 


0.13 


0.003 


62 1553 + 24 


35.28 


1 






UGC1115 


30.57 


1.03 






O2208 


40.43 


16.67 


X 


X 



Table 3. IRAM 30m fluxes at 3 mm and 1 mm. X means that 
there was no data, and " ..." refers to undetected flux. 



ing regions and/or an AGN, and the cooler components (A > 
100/im) associated with more extended d ust heated by the in- 
terstellar radiation field (iKennicutj 119981) . We complemented 
the 30m data with the IRAS fluxes taken from the Nasa Extra- 
galactic Database (NED) for the 60/mi and 100/mi wavelengths. 
From this sample, a total of 38 galaxies were detected by IRAS. 
18 of them are considered a clear detection; 15 galaxies are up- 
per limits for both 60/iin and 100/im wavelength and 5 galaxies 
are upper limit only at 100//m and detected at 60//m. For that 
study we will use hereafter the 18 galaxies that were clearly 
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detected in both wavelengths, the galaxies and their fluxes are 
listed on Table [4] 



VrdldXy IX dlllc 


fw /^m 


/ioo A/m 




(nuy) 




(mjy) 






4jJ ± 


A^ 1 

CD. j 


1 A7^ 


+ Zjl 


oo 


loU ± 


L 1 


C 1 A _i_ 


Izz 


Ill 1 
jL-iv Z /Z. 1 


JJO + 


S3 A 


i no a 


± 1j4 


11 A 


j4o ± 


oi.y 


£ £0 _i_ 


oq o 
OJ.y 


JKJJ 


ZVo + 


AA H 
44. / 


4jU ± 


D /. J 




lUO/ 


j_ 1 ao 


vol ± 


\AA 

144 


'ill 


o /u ± 


1U1 


J 1 1 ± 


JJ. / 


jLK 4UZ 


257 ± 


38.6 


1052 


± 158 


NGC 315 


368 ± 


55.2 


460 ± 


69 


NGC 708 


200 ± 


34 


660 ± 


157 


NGC 4278 


580 ± 


53 


1860 


±60 


NGC 6251 


188 ± 


28.2 


600 ± 


90 


NGC 7052 


524 + 


78.6 


1150 


± 173 


B2 0116+31 


150 ± 


22.5 


524 ± 


78.6 


B2 0648+27 


2633 


±395 


1529 


±229 


B2 0836+29B 


472 ± 


70.8 


595 ± 


89.2 


B2 1101+38 


181 ± 


22 


361 ± 


68 


OQ 208 


753 ± 


113 


1029 


± 154 



Table 4. IRAS data for the galaxies that have been detected at 
60 yum and at 100 yum. 



3. Molecular Gas 

3.1. Molecular gas mass 

In Table [5] we present the results of the observations, where 
Ico(i-D) and Ico(i-i) are the integrated intensities (J TmbcIv) 
for the 12 CO(1-0) and 12 CO(2-l) line respectively, 5I C o(\-Q) 
and 5Ico(2-\) are the standard error on /co(i-O) and Ico(2-\), 
V 'width coo-0) and V width C o{2-\) are the velocity width of CO(l- 
0) and CO(2-l) respectively; and finally, M# 2 is the calculated 
molecular gas mass for e a ch g alaxy in the sample. 
We follow iGordon et alJ <[l992) to derive Hi masses from the 
12 CO(1-0) line observations. Our temperature unit is expressed 
in T* A antenna temperature scale, which is corrected for atmo- 
spheric attenuation and rear side-lobes. The radiation tempera- 
ture Tr of the extragalactic source is then: 



4 A j K Ta 
n D ?] A Ll s 



(1) 



where A is the observed wavelength (2.6 mm), D is the IRAM 
radio telescope diameter (30 m), K is the correction factor for 
the coupling of the sources with the beam, t/a is the aperture 
efficiency (0.55) at 1 15 GHz, Ta is an antenna temperature which 
is F e ffT* A in the IRAM convention, explicitly T A = Q.92T* A , and 
L~ls is the source solid angle. 

Without taking into account cosmological correction, because 
of the low redshift and using the standard CO-to- H? conversion 
factor, 2.3 x \Q 10 cmr 2 {K km s~ l )~ x suggested by Stron getal] 
( 1988), the column density of molecular hydrogen can be written 

as: 



N(H 2 ) = 2.3 x 10 



2d 



f 

J I in 



T^dv{mol.cm 2 ) 



(2) 



Where dv is the velocity interval. In Eq. \T\K = Q,/Qs is the 
factor which corrects the measured antenna temperature for the 
weighting of the source distribution by the large antenna beam in 



case of a smaller source. We have defined the source solid angle 
as: 



-J 

*J SOI 



0(0, \p)dL~l 



(3) 



where 4>(6, if/) is the normalized source brightness distribution 
function. The beam-weighted source solid angle is 



*J SOI 



(4) 



where f(Q, <A ) denotes the normalized antenna power pattern 
(lBaarslll973h . Experiments have shown that we can approxi- 
mate / by gaussian beam. An exponential law of scale length 
h = Dg/10 is taken to model the source distribution function. 
The factor 10 is used to estimate the scale length of the molecu- 
lar gas distribution from the optical diameter (B band). We used 
our own data from the CO(1-0) map of 3CR 31 and the CO(1-0) 
map of elliptical galaxies presented by Young (2002). In EqO 
that optical diameter is represented by # s .Usuing the asumption 
that the gas surface density is fi(r) oc e^^ h and as long as the 
source is smaller then the beam size, we have the following: 



J~9J2 
sin(e)e~<>s d6 

sin{ff)e " s V 



K = 



Jo 



dO 



(5) 



If Ico is the velocity-integrated temperature for the I2 CO(1-0) 
line in T* A scale and given the IRAM-30m parameters, the total 
mass of H2 is then given by 

M Hl = 5.86 x 10 4 D 2 KI co (Mb) (6) 
with the distance, D, in Mpc. 

Using the survival analysis statistics software (ASURV), which 
takes into account the upper limits, besides the detections, we 
calculated a median value for the molecular gas mass of 2.2 X 
10 8 M (and an average of 3xl0 8 M ) for our complete sample 
of 52 galaxies. Their molecular gas mass distribution can be seen 
on FigQ] The molecular gas distribution is normalized to one for 
the whole sample including the detections and the upper-limits 
for CO(1-0) observations. On the same Figure we show in darker 
color the histogram of only the galaxies detected, the median 
value of the molecular gas mass of these 30 detected galaxies is 
then 3.4 x 10 8 M . 

To normalize the sample we divided the distribution per bin of 
the histogram by the total number of galaxies in the sample, and 
to normalize the subsample, that includes the detected galaxies 
only, we divided the distribution per bin by the total number of 
galaxies in the complete sample, with the purpose of comparing 
both samples. The range of the molecular gas mass (in Log) for 
the sub-group of the detected galaxies is between 6.3 and 10.2. 
The number of galaxies with a molecular gas mass lower than 
6.3 (in Log) is obviously unknown. 

The molecular gas mass in the sample ranges over 4 orders of 
magnitude, from 3 x 10 6 M for 3C272.1 to 1.4 x 10 10 M Q for 
OQ 208. Most of the galaxies have a molecular gas mass be- 
tween 8.3 (2 x 10 8 M ) and 8.7 (5 x 10 8 M ), this is true for the 
complete sample including the upper limit calculation as for the 
subsample that includes only the detected galaxies. 
If we look at the relation between the optical luminosity (L^) 
and the molecular gas mass of the galaxies (M//,), as shown in 
Fig. |2 we find that there is no correlation, which is interpreted 
as that the molecular gas mass is independent of the size of the 
galaxies. The FIR Luminosity (which is more detailed in sec- 
tion |5}, and the molecular gas mass, normalized by the optical 
luminosity has also no correlation, as can be confirmed on Fig. 
® 
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Table 5. Molecular Gas data 
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The * means that this galaxy is a tentative detection, the ** means this galaxy is detected. The letter a is for (l-0)&(2-l) detected galaxies, b is for 
the CO(l-O) only detected galaxy and c is for the galaxies detected in CO(2-l) only; note that in this cases the velocity with is represented as — 
since the CO(l-O) value was calculated using the average of the ratio and therefore the width of the line is unknown. The case of NGC4278, that 
has an a* is because this galaxy has been detected in (2-1) and tentatively detected in (1-0). In M Hl the * have the same meaning, only that in this 
case refers strictly to CO(l-O) since the integrated velocity for 3 mm was the one used to calculate the molecular mass. 



3.2. CO(2-1)-to-CO(1-0) line ratios 

The line ratio between the CO(2-l) and the CO(l-O) transitions 
is computed by comparing the integrated intensity ratio of the 
lines Ico(2^i) / Ico(i^O) where the intensity was measured on 



one point at the center of each galaxy. Because of the different 
beam sizes at CO(l-O) and CO(2-l), the intensities should be 
compared only after having carried out small maps in CO(2-l) 
to sample the CO(l-O) beam. However, in most of our galaxies, 
the CO emission is expected to be almost a point source, at both 
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Fig. 1. This figure represents the histogram of the sample. 
The lighter color represents all the galaxies in the sample. Using 
ASURV we were able to calculate the mass of the sample using 
both, upper limits and detected galaxies. Plotted at the bottom, 
with the darker color, are the detected galaxies only. 









• •• . ' 




• 




















• 












• • 












• * • . 




• 








• 

• 




• • • 








• 


• 










• 







10.5 



Fig. 2. Molecular gas mass versus blue luminosity with the upper 
limits for the molecular gas mass indicated. 



frequencies. Except for NGC 4278, the distance of our detected 
galaxies is around or above 100 Mpc, where the CO(l-O) 
beam is 10 kpc, and in elli ptical galaxie s , the g aseous disks are 
expected to be, according to Youn geTal] ([2008), concentrated in 
the nuclear disks at Kpc scales. Therefore, equal Ico intensities 
in both lines will result in ratios of ~ 4, the point source dilution 
ratio of the two beams (to have a more detailed derivation of 
beam couplings, see AppendixlBl). 

As previously noticed by Lim et a D (EOOO), the 2 galaxies 
studied in their paper have a stronger observed intensity in 
CO(2-l) than in CO(1-0). This sample has a line ratio well 
over unity. Seven of our galaxies have been detected in both 
frequencies and 6 more have been detected only in CO(2-l). 
From the galaxies detected in both frequencies, only NGC 
7052 has an integrated intensity stronger in the CO(1-0) line 
than in the CO(2-l) line. The maximum line ratio was found 
to be 3.1 (for NGC 4278) and the average value is 2.3 + 0.1. 




-9.8 -9.6 -9.4 
LoglLfir/Lb) 



Fig. 3. Molecular gas mass versus FIR luminosity, both normal- 
ized by the blue Luminosity, using only detected values. 



Galaxy Name 


Ratio 


3CR31 


2.5 ± 0.2 


3CR 264 


2.1 ±0.2 


3CR 272.1 


2.5 + 0.4 


NGC 326 


2.8 ± 0.3 


NGC 4278 


3.1 ±0.4 


NGC 7052 


0.8 ±0.1 


B2 0116+31 


2.4 + 0.3 


B2 0648+27 


2.0 ± 0.3 



Table 6. CO(2-l)-to CO(1-0) line ratio of the detected galaxies, 
not corrected by the different beam dilutions. 



When correcting by the factor 4 of the beam dilution ratio, the 
average value is ~ 0.6. This corresponds to an average over 
a moderate density disk, where most of the CO emission is 
optically thick, but sub-thermally excited. The maximum ratio 
obtained precisely for the most nearby galaxy NGC 4278 is 0.8, 
comparable to the t ypical value obta i ned fo r the central parts of 
spiral galaxies bv lBraine & Combesl (Il992h . 
Figure|4]is a plot of the integrated intensity Ico(2^\) vs. /co(i^O) 
with a line fit indicating a clear correlation. The dotted line 
in the plot is for the intensities with the same value for both 
lines. Table [6] lists the line ratios for each galaxy and to see the 
spectra of those galaxies detected in both frequencies see the 
online material, where the spectra of all galaxies in the sample 
are plotted. Note that the galaxy B2 0116+31 presents in its 
CO(1-0) line profile a strong absorption line, as well as the 
double line profile that can still be clearly seen. This absorption 
is the signature of molecular gas mass on the line of sight 
towards the AGN covering a very small area. Most of the radio 
galaxies observed have a radio continuum in millimeter strong 
enough to be detected in absorption, but only in one galaxy the 
CO(1-0) transition is absorbed towards the radio continuum. 
The double horn profile of B2 01 16+31 is visible as well in the 
CO(2-l) transition line, where the absorption is much weaker, 
as well as the continuum at this frequency. The absorption line 
in this galaxy should cause an underestimation of the CO(1-0) 
integrated intensity and therefore an overestimate of the line 
ratio. 

For those galaxies detected in CO(2-l) emission line, and 
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3.3. Fanaroffand Riley classification 



1.0 




Linear Fit 
Ico(l-O)=lco{2-1) 



0.0 

Log lco(l-O) 



Fig. 4. Plot of the line ratio of the galaxies. The point that is 
outside of the linear fit is of the galaxy NGC 7052 



not detected in CO(l-O) we derived a value for the integrated 
velocity /co(i-o) using the ratio of 2.3 calculated with the 
galaxies detected in both lines. Again, for a closer view to the 
spectra of those galaxies please see the online material. The 
values of the derived molecular gas masses are already included 
on Table [5] 



3.2.1 . Origin of the low line ratio 

As mentioned previously, the observations were done only in 
one point, at the center of the galaxies. However, the sources in 
the present sample are distant, and the expected extent of the 
CO emission is almost that of a point source, with respect to the 
large CO beams. 

When we can appreciate from the HST optical images or from 
the CO interferometric map of 3CR 31 (see Fig. \C.2\ . the size 
of the molecular gas emission, it is possible to correct for the 
beam/source coupling ratio. 

Taking into account the only 2 galaxies from this subsample for 
which we have the size of the molecular gas emission - 3CR 
31 (§E4> and NGC 7052 (§|A), we could roughly say that the 
size of the galaxies varies from 4" to 8", so that the correction 
factor for this subsample would be between 0.25 and 0.35. For 
3CR 3 1 the result is more accurate since we have the CO map 
from the PdBI, and as can be derived from the map, we can 
see the galaxy has a molecular gas size of about 8", implying a 
correction factor of about 0.3. If the line ratio of 3C R31 is 2.5, 
then after the correction, the ratio would be 0.8. This value is 
lower than one and th erefore it could ind i cate th at the CO line is 
sub-thermally excited Br aine & Combes! (1 19921) 
In average, we find a value of ~ 0.6 for the CO(2-l)-to-CO(l-0) 
line ratio. 

This is compatible with a disk where the CO line is optically 
thick and sub-thermally excited, as is the disk of spiral galaxies 
in general. Only when it is possible to resolve the center of the 
disk, in nearby galaxies, the ratio has been observed to rise up 
to 0.8-1 dBraine & Combesll 1992b . 



Radio galaxies are divided mainly in two groups, the Fanaroff 
and Riley type I (FR-I) and type II (FR-II). As explained by 
Fanaroff & Rilevl d!974l) . the sources were classified using the 
ratio of the distance between the regions of highest brightness of 
the radio continuum on opposite sides of the central galaxy or 
quasar, to the total extent of the source measured from the low- 
est contour; those sources from which the ratio was less than 0.5 
were placed in class I, and those for which the ratio was greater 
than 0.5 were placed in class II. There is a third classification 
called FR-c galaxies. These galaxies have a very compact radio 
continuum emission. Their radio morphologies suggest that they 
are compact versions of the classical FR-II's, although why they 
are so small is not yet established: It is hypothesized that these 
are ei ther young FR-IFs or FR-IFs trapped in a dense environ - 
ment (Fan ti et al .111 990t [C? Pea et al ] [l99lllFanti & Fantilfl99i . 
In the present sample, 69.2% are FR-I type galaxies; 19.2% are 
FR-II type galaxies and 1 1.5% are FR-c type galaxies. 
We have previously said that the molecular gas mass average of 
our sample is 2.2 x 10 8 M . FR-I and FR-II type galaxies do be- 
have differently from one another. The FR-II type galaxies have 
the largest molecular gas mass compared to the FR-I type galax- 
ies and the FR-c type galaxies, as shown in Table [7j using the 
detected galaxies and the upper limits with the survival analysis 
statistics. The mean molecular gas mass is lower for the FR-I 
type galaxies than for the FR-II type galaxies with 1.7 x 10 8 M 
and 8.1 x 10 8 M respectively. This difference is clearly visible 
in Figure [5] where the molecular gas mass distribution is shown 
according to the types of radio galaxies: FR-I, FR-II and FR-c. 
We can draw the following conclusions regarding the different 
types of radio galaxies and their molecular gas masses: 

1 . For the FR-I types, which are the less powerful AGN, the el- 
liptical host galaxy does not need much molecular gas mass 
to host the radio AGN. Molecular gas masses can be as low 
as 10 6 M , which is a few Giant Molecular Clouds (GMC). 

2. For the median value of the molecular gas mass, FR-II type 
galaxies (8.1 x 10 8 M ) are clearly more massive than the 
FR-I (1 .98 x 10 8 M ) and FR-c (2.0 x 10 7 M ) type galaxies. 
Note that there are only a few FR-c galaxies and therefore 
the statistical values for them are not very reliable. 



Evan s et al ] (120051) found only an 8% detection rate of CO 
emission in their FR-II galaxies, a very low value compared to 
the ~ 35% detection rate of the FR-I detected in their sample. 
We also found a higher detection rate in FR-I galaxies compared 
to the FR-II galaxies. In our sample, the detection rate is 38.5% 
for FR-I, 11.5% for FR-II and 7.7% for FR-c galaxies. Note 
that the FR-II being more powerful AGN, they are rarer and are 
found at larger distances. This can also explain both their lower 
detection rate, and their higher median molecular content. 
To better estimate the Malmquist bias in our sample, we plot in 
Figure |6]the molecular gas mass vs. the redshift (z). It is clear 
that for higher z the galaxies tend to have higher molecular gas 
mass, following the sensitivity of the telescope. The sensitivity 
limit was computed, assuming a typical value of 300 km/s for 
the velocity width and 1 mK for the rms noise temperature. 
It is clear that there is a larger number of FR-II galaxies, 
compared to FR-I and FR-c galaxies, at higher z, implying a 
larger threshold of the upper limit for the FR-II type galaxies. 
This also agrees with the idea that FR-II galaxies are stronger 
AGN and more luminous, and this is also why they are seen 
at a higher distance. A factor 3 in the median mass can be 
expected, and we conclude that the difference between the radio 
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I FR-c 
I FR-II 
I FR-I 



Galaxy 






(XlO M ) 


3CR 31 


16.81 ± 2.38 


3CR 264 


3.37 ± 0.74 




A 7^ _i_ 1 1 Q 


3CR 449 


J. J7 IT- !_/. / J 


NGC 3801 


5.49 ± 1.63 


NGC 7052 


1.96 ± 0.41 


B20116+31 


60.63 ± 16.92 


OQ 208 


136.51+ 17.94 


Average Value: 


29.14 ±5.24 



nHMI 



Fig. 5. Molecular gas mass distribution in the radio galaxies de- 
pending on their Fanaroff-Riley classification. 




• »FR-I 
A A FR-I I 
X X FR-C 
Theoretical 



Lower Limit 



0.06 

z 



Fig. 6. M//, versus z and a lower limit of the theoretical value for 
the mass with respect to the distance. 



galaxy types for the molecular gas mass could be only due to 
the Malmquist bias. 



Table 8. Molecular mass gas of the galaxies that present the dou- 
ble horn profile. 



We also note that 4 of our 8 galaxies have been detected in both 
transitions (3CR 31, 3CR 264, NGC 7052 and B2 1 1 6+3 1 ) and 
in all 4 cases the CO(2-l) line is stronger than the CO(1-0) line 
(more details about the line ratios of these galaxies in § 13.21 ). 
The galaxies that present the double horn profile feature are 
also within the largest H2 masses. The average value for this 
subsample is 29.1 x 10 8 M , an order of magnitude higher value 
than the average value of the complete sample for the detected 
galaxies only (2.2 x 10 s M ). 

We note that in some cases the absorption of the CO transitions 
could mimic a double horn profile. It appears from the peak 
temperature in the double horn profiles that in all cases, except 
for the CO(2-l) emission in 3CR 31, the double horn profile 
could be perturbed by an absorption line. 



4. Dust 

Dust is an important ingredient of the ISM, associated with the 
molecular gas. It is heated by stars and/or an AGN, and its FIR 
emission depends strongly on its temperature, giving informa- 
tion about star formation (SF) in the galaxies. We analyze the 
012 dust emission at 60 //m and 100 yum in our sample. The fluxes 
at these frequencies give an estimation of the temperature and 
the mass of the warm dust component. On a companion paper 
we will use the fluxes at 12 /mi and 25/vm as well as the fluxes 
from the Spitzer satellite using both cameras, IRAC (3.6, 4.5, 
5.8, 8 11) and MIPS (24,70, 160 //) to create the Spectral Energy 
Distribution (SEDs), and to study in more details the synchrotron 
and thermal emission in each galaxy of the sample. 
The dust temperature T d is calculated by adopting a modified 
black-body emissivity law: 



3.4. Molecular gas disk 

The double horn CO line profiles observed are characteristic 
of an inclined molecular gas rotating disk. In the present set 
of data, the double horn feature is present in 8 galaxies: 3CR 
31, 3CR 264, 3CR 403, 3CR 449, NGC 3801, NGC 7052, B2 
0116+31 and OQ 208, i.e. 29% of the detected galaxies in our 
sample. In the online material you can find the spectra of those 
galaxies that present a double horn profile. In the case of 3CR 
3 1 and 3CR 264 the double horn profile is clearly visible in both 
transitions, CO(1-0) and CO(2-l), and for 3CR 403 and OQ 208 
we see only this feature in the CO(1-0) transition since there 
is no data for the CO(2-l) transition. Five of the galaxies that 
present the double horn profile are FR-I types. Only one galaxy, 
3CR 403, is an FR-II type, and two galaxies, OQ 208 and B2 
0116+31, are FR-c type. 



Sy OC KyBy(T d ) 



(7) 



k v is the dust emissivity and B v (T d ) is the Planck function at the 
dust temperature T d . Then we calculate the dust mass as follows: 



M dust = 4.78/ 10 o£>£(e (143 - 88/rrf) - 1) M e 



(8) 



Where /ioo is the flux at 100 fim in Jy, Di the luminosity distance 
in Mpc and T d the dust temperature in K . The median value of 
the dust temperature for this sample, using the galaxies detected 
by IRAS, is 35. 8K and the median value of the dust mass for the 
sample, only for the galaxies detected by IRAS, is 2.5 x 10 M . 
The detailed list of values for each detected galaxies at the fre- 
quencies of 60 yum and 100 fim is given in the Table|9] 
The Mh 2 /Mdust gas-to-dust mass ratio for this sample, using the 
galaxies detected with both IRAS and IRAM-30m, is 260. On 
Figure [7] the distribution of the gas-to-dust mass ratio is shown, 
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Galaxy Name 


Detected & Upp 
Mean 


er Limits 
Median 


Detected ( 
Mean 


)nly 
Median 


FR-I 
FR-II 
FR-c 


r 7 v rn 8 +2.3xio 5 

l,XlV _i.2xl08 
8 1 v lf> 8 +1.22x10' 

2.1 X 10 s +9 - 2x10 * 

-4.8x10? 


1.86 x 10 8 
4.45 x 10 8 
2.03 x 10 7 


2.3 x 10 8 

-1.6x10 s 
Llxl0 9 +1.7X1* 
-6.5x10 s 

8.5 x io 8 + ;f ;° 


2.01 x 10 8 
4.6 x 10 8 
1.2 x 10 9 


Total 


? 3 X 10 8 +3-2X10 5 


2.2 x 10 8 


1 1 y in 8 +5.2x10" 


3.4 x 10 8 



Table 7. Mean and median values of the molecular gas mass of all subsamples separately in logarithmic values in units of M ( 



4.C 4.5 5.0 5.5 6.0 6.5 



7.0 7.5 



0.45 r 
0.40 - 
0.35 
0.30 - 
& 0.25 - 
1 0.20 - 
0.15 
0.10 
0.05 - 



0.01 



°>0 



T,,,.,V<1 



Fig. 8. Histo; 
temperature. 



a) b) 

;rams of the dust characteristic of the sample. Figure a) is the dust mass in Logarithmic values and Figure b) is the dust 
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3.5 - 

3.0 - 
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I 2.0 - 
jj 

1.5 - 
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0.5 - 
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Fig. 7. Histogram of the Gas-to-Dust mass ratio. 

where the galaxies range between 20 and 4450 in the gas-to-dust 
mass ratio. Comparing with the sample that includes the upper 
limits for the molecular gas mass gives us an upper limits for 
the ratio of 254 on average. Figure [8^ is the histogram of the 
dust mass of this sample. Note that most of the galaxies in this 
sample have a dust mass in the range 7.4 - 30 x 10 5 M with 
80% of the detected galaxies having a dust mass greater than 
7.4 X 10 5 . The dust temperature is mostly between 27.0 and 34.4 
K (see Figure [8] b) although there are 2 galaxies with an esti- 
mated dust temperature of 64 and 71.4 K respectively. These 2 
galaxies with a hotter dust component than the others are 3CR 
327 and B2 0648+29B. These temperatures are derived assum- 
ing that the heated dust is mainly radiating by the FIR traced by 
the IRAS fluxes at 60 and 100 fim . Nevertheless we cannot dis- 



Table 9. Dust mass and temperature. 



Galaxy Name 


Log(M dusl ) 


Tdust 




(M ) 


(K) 


3C31 


6.8 


28.3 


3C88 


7.1 


27.0 


3C272.1 


4.6 


36.4 


3C274 (M87) 


4.1 


47.4 


3C305 


6.4 


39.5 


3C321 


7.1 


50.8 


3C327 


6.4 


71.3 


3C402 


7.0 


27.8 


NGC 315 


5.5 


42.9 


NGC 708 


6.3 


29.7 


NGC 4278 


5 


20 


NGC 6251 


6.9 


30 


NGC 7052 


6.2 


34.1 


B20116+31 


7.4 


29.2 


B2 0648+27 


6.2 


68.5 


B2 0836+29B 


6.8 


42.8 


B2 1101+38 


6.2 


35.3 


OQ 208 


7.2 


41.2 



card a contribution of warmer and colder dust components than 
traced by IRAS. This issue will be addressed in a companion 
paper. 

5. FIR vs. CO 

5.1. FIR & CO Luminosity 

We calculated the CO(1-0) line luminosity, L' co , which is ex- 
pressed as velocity-integrated source brightness temperature, 
T/,AV, times source area Q. S D^, and which for our sample re- 
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suited in a mean value of 4.2 x 10 7 K kms~ 1 pc 2 . This value was Table 10. Values of CO(l-O) Luminosities and FIR Luminosities 
calculated using ASURV for the upper limits. The CO(l-O) lumi- 
nosity mean value is 6.4 xlO 7 K kms~ l pc 2 calculated only with 
the detected values of the CO(l-O) e mission lines. The for mula 
used for the calculation is taken from Solo mon et a fl <fl997h : 



L' co = 23.5 xn 2 stb Dll C o(l+ zT 



(9) 



Where Q. s *t, is the solid angle of the CO(l-O) emission for the 
source convolved with the telescope beam, Di the distance in 
Mpc, Ico is the integrated velocity in K km/s (in T„,h) and z the 
redshift of the galaxy. If the source is smaller than the beam, 
which is the case for all the galaxies in this sample, we can as- 
sume that Q. s *t> ~ Qb which is given in arc sec 2 . 
Lfir is the far-IR luminosity using the fluxes at 60 pm, associ- 
ated with warm component, and 100 pm. associated with cooler 
components. It was com puted using the following relation from 
San ders & Mirabelld 19961) : 



^FIR 



1 + 



ft 



11)0 



2.58/ e 



60 



^60 



(10) 



Where fa) and /ioo are the fluxes at 60 and 100 pm respectively 
and are in units of Jy; L(,q is the luminosity at 60 /im in L Q and it 
is represented as: 



Log(L 60 ) = 6.014 + 2Log(D) + Log(f 6Q ) 



(11) 



Where D is the distance of the galaxy in Mpc. 
The results can be seen in Table [TOl for which the median of the 
Lfir = 6.9 x 10 9 L Q . This value was calculated using ASURV 
to take into account the upper limits. 

5.2. Star Formation 

The ratio of Lfir/L' co is normally related in spiral galaxies to 
the star formation efficiency (SFE). This assumes that most of 
the FIR emission is coming from dust heated by young stars. 
In the case of powerful AGN, there is always the possibility of 
an added FIR contribution from dust heated by the AGN in a 
dusty toru s. However, this co ntribution is rarely dominant (see 



i sty tor 

eg. lGenzel & C esarskv 2000). Figure [9] is a plot of the galaxies 
in our sample that represents the Lfir vs L' co for the detected 
galaxies in both FIR and CO with their error bars, the galaxies 
with upper limits, represented with arrows (eg, down arrow is for 
the CO upper limit and FIR detection, right arrow is for the CO 
detected and upper limits for the FIR and the diagonal arrows 
are for the CO and FIR upper limits) the green line with the two 
circles are for when the Lfir was calculated with a 60/mi flux de- 
tection and 100/mi flux upper limit, representing the lower and 
upper limit for the Lfir. 

The linear relation between Lfir and L' co is log(L' co ) = 
0.8log(L F m) + 0.1 with a x 2 of 0.024. Finally the relationship 
itself has a median value for the Lfir/L' co ratio of about 52 
L e (K km s- 1 pc 2 )- 1 . 

We note that the detection rate for both FIR and CO together 
is independent of the type of radio galaxy. Out of the 15 radio 
galaxies detected for both FIR and CO, 8 (22%) are FR-I (3CR 
31, 3CR 272.1, 3CR 274, 3CR 305, NGC 315, NGC 708, NGC 
7052 and B2 0836+29B), 3 (30%) are FR-II (3CR 88, 3CR 321 
and 3CR 327) and 4 (66%) are FR-c (NGC 4278, B2 01 16+31, 
B2 0648+27 and OQ 208). 

We note as well that from the 15 galaxies detected for both FIR 
and CO, only 4 have a double horn profile feature (3CR 3 1 , NGC 
7052, B2 0116+31 and OQ 208). 



Galaxy 


L °g( L 'co) Log(L Fm ) 
(K km s- 1 pc 2 ) (L ) 


3C31 


8.51 


9.77 


3C40 


<7.86 


<9.30 


3C66B 


7.23 


<9.32 


3C83.1 


7.78 


<9.50 


3C88 


7.62 


9.94 


3C98 


<8.12 


<9.54 


3C129 


7.38 




3C264 


7.81 


<9.46 


3C270 


<7.13 


< 8.51 


3C272.1 


5.77 


8.23 


3C274 


7.01 


8.25 


3C296 


<7.99 


<9.53 


3C305 


8.60 


10.21 


3C321 


9.12 


11.43 


3C236 


<9.10 


< 10.28 


3C327 


8.69 


11.24 


3C353 


7.80 




3C386 


7.52 


< 8.75 


3C402 


<8.07 


9.90 


3C403 


8.43 


< 10.62 


3C433 


<9.17 


< 11.20 


3C442 


7.08 




3C449 


7.84 


<9.34 


3C465 


<8.14 


<9.66 


NGC 315 


7.20 


9.51 


NGC 326 


7.98 


<9.88 


NGC 541 


7.64 


<8.79 


NGC 708 


8.01 


9.38 


NGC 2484 


< 8.54 




NGC 2892 


<8.10 




NGC 3801 


8.03 


<9.60 


NGC 4278 


6.00 


8.08 


NGC 5127 


7.14 


<8.45 


NGC 5141 


< 7.97 


<8.72 


NGC 5490 


< 8.00 


<8.74 


NGC 6251 


< 8.09 


9.69 


NGC 7052 


7.58 


9.65 


B2 0034+25 


<8.27 




B2 0116+31 


9.06 


10.39 


B2 0648+27 


8.37 


11.03 


B2 0836+29B 


8.96 


10.76 


B2 0915+32B 


<9.10 


<9.93 


B2 0924+30 


7.59 




B2 1101+38 


< 8.37 


9.75 


B2 1347+28 


8.50 




B3 1357+28 


< 8.90 




B3 1447+27 


< 8.16 




B3 1512+30 


<9.23 




B3 1525+29 


<9.06 




B3 1553+24 


< 8.72 




UGC7115 


< 8.19 




OQ 208 


9.41 


11.12 



As mention, Lfir/L' co is interpreted as a measure of SFE. Then, 
the Star Formation Rat e (SFR) per unit gas mass, according to 
iGao & Solomon! d2004l) . is represented as: 



Msfr * 2 x 10- 10 {L FI RlL )M &yr 



(12) 



From this estimation we compute a mean SFR for this sample 
of about 3.7 ±1.5 M e /yr, assuming that all the FIR emission 
comes from the dust heated by the young formed stars. This is 
likely an overestimation since it does not take into account the 
AGN heating of the dust. 
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L 7 H 111 I j 1 \_- 






J null 1 ratio 
L-'FIRI '-'CO iaLlu * 




(Ad) 


(Kkm/spc 2 ) 


(L Q (Kkm/spc 2 )- 1 ) 


Xhis sample 


8.9 x 10 y 


4.2 x 10' 


52.0 


Wiklindetal. (1995) 


6.5 x 10" 


1.6 x 10 s 


39.4 


Mazzarella et al. (1993) 


1.94 x 10" 


1.1 x 10" 


82.2 


Evans et al. (20051 


3.6 x 10 10 


1.2 x 10" 


35.4 


Bertram et al. (2007) 


4.2 x 10 10 


9.2 x 10 8 


49.1 


Solomon et al. (1997) 


1.2 x 10 12 


7.9 x 10" 


146.8 




Detected 
Upper Limits 

flux60 Detected, F100 UL 

— Linear Fit 

- - Linear Fit - all samples 



5.0 8.5 9.0 



9.5 10.0 



Fig. 9. Lfm vs. L' co for the galaxies in this sample. The upper 
limits are shown, where the Leo upper limits are represented by a 
down arrow, the Lpm upper limits are represented by a lef arrow 
and the Lco-Lfir upper limits are represented by south weast 
arrows. It is also possible to appreciate the galaxies in which the 
flux at 100 /im was an upper limit and the flux at 60 pun was 
detected eith their error bars, in case Leo was detected and the 
down arrow in case Leo was upper limits 



5.3. CO-FIR relationship 

In systems with ongoing star formation (SF), the light from both 
newly formed and ol der s tars can be absorbed by dust and repro- 
cessed into the FIR. iBelll (2003) states that since young stars in 
HII regions heat up dust at relatively high temperatures (with a 
low 100 jum-to-60 yum ratio of about ~1) we could have an idea 
of where is this IR emission coming from, and this ratio leads 
to a wide range in/ioo / j m //6o J um on galaxy- wide scales, from ~10 
for early type spiral galaxies to < 1 for the most intensely star 
forming galaxies suggesting that earlier types are influenced by 
old stellar population. 

Our sample has a median value of /loo/jm/Zeo/jm = 1-9 with its 
smaller value being 0.55 for 3C 327, which could be interpreted 
as galaxy with high temperature heated up by young stars in HII 
regions. And its highest value is 4.5 for 3CR 88, which is, still, a 
small value c ompared w ith ~10 for quiescent early type spirals 
suggested bv lBelll (120031) . 

Fig. |9]shows a linear relation in our data for the L' co and Lpm. It 
is evident that the CO luminosities increase linearly with increas- 
ing Lfm with a relation such as log{L' co ) = QMogiLpm) - 0.6. 
On the plot we see there are two linear fits, the first one, the 
solid line, represents the linear fit of the galaxies from this sam- 



I 1.5 



A AFR-I 

FR-II 
O OFR-C 



o 



o 



A 
A 



10.5 



Fig. 10. Star formation efficiency compared to the blue luminos- 
ity of this sample of galaxies. 



pie only. On section [6] we will compare our sample with other 
samples, this second fit, the dashed line, is the fit that we got 
from fitting all the galaxies together. 

We have related the line ratio of the galaxies clearly detected 
in both lines with the gas-to-dust mass ratio as can be seen in 
the online section on Figure IB. II From this we see that most of 
the galaxies have a gas-to-dust mass ratio lower than the median 
value for the ratio of the complete sample. Recall that the aver- 
age gas-to-dust mass ratio of the sample is ~300. 
The galaxy 3CR 88 has the lowest gas-to-dust mass ratio, 17, a 
very low value compared to NGC 4278 with a ratio of 600. 
We also compared the SFE (L'/Lfir) vs. Lb, as shown in figure 
[10]where we found no correlation at all. In the figures the trian- 
gles are for the FR-I type galaxies, the squares are for the FR-II 
and the circles for the FR-c. The sample all together do not show 
correlations as well as all the subsamples separately. 



6. Comparison with other samples 

6. 1 . General description of the comparison samples 

We compare our sample with other samples selected under dif- 
ferent criteria. 



- lEvans et all (120051) observed a sample of elliptical radio 
galaxies, like our own, but selected with the IRAS fluxes 
densities at 60 //m, fbo or 100 p/m, /ioo^m, greater than 
0.3 Jy, since they assumed that this level is a good indicator 
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of a dusty, gas-rich, interstellar medium. 

- IWiklind etal] d 19951) observed a sample of elliptical galax- 
ies selected on the basis of their morphological type. They 
argue that there is a major difficulty in defining a sample of 
elliptical galaxies due to uncertain morphological classifi- 
cation. According to them, the best criterion for defining 
an elliptical galaxy appears to be the Vaucouleurs r 1 ^ 4 
luminosity profile, and based on this, they selected genuine 
ellipticals, or galaxies that had a consistent classification as 
E in several catalogs. A 100 /vm flux > 1 Jy and 5 > 0° were 
an additional criterion. 

- iMazzarella et alj d 19931) studied an infrared limited sample 
with 60 /mi flux greater than 0.3 Jy. The objects were chosen 
to fully cover the sky, and to be relatively near (z < 0.1), 
with no other selection criteria. 

- iBertram et al.l (120071) studied a sample of galaxies hosting 
low-luminosity quasi-stellar objects (QSOs) where they 
argue that an abundant supply of gas is necessary to fuel 
both, an active galactic nucleus and any circumnuclear 
starburst activity. The only criterion of selection was their 
small cosmological distance: only objects with a redshift 
z<0.060 were chosen. 

- ISolomon et al.l d 19971) studied a sample of 37 infrared- 
luminous galaxies in the redshift range z=0.02-0.27. Eleven 
of these galaxies have a 60 /mi flux $ 60 > 5.0 Jy, and are 
part o f the near bright galaxy sample dSanders et al.lll988l 
1 199 lh. Twenty gala xies were chosen from a redshift survey 
dStrauss et alJll992h of all IRAS sources with 60 //m fluxes 

> 1.9 Jy. And a few lower luminosity sources were also 
included in the sample. 



6.2. Molecular gas mass 

Comparing our sample with the other samples mentioned above, 
we noticed that Wikli nd et all dl995l) and our sample have an 
average of the molecular gas mass of a few 10 8 M (using 
only detected galaxies, not the survival analysis median value) 
and that both samples a re being hos t ed by ellip tical galaxies; 
Mazz arella et al.l (Il993l) . lEvans etal.1 (120051) and IBertram et all 
(2007) samples all have a median value of the molecular gas 
mass about 10 9 M , their samples are FIR selected galaxies 
or galaxies in interaction, therefore with an expected higher 
molecular gas mass than in the sample hosted by normal 
elliptical galaxies. They conclude that the majority of luminous, 
low-redshift QSOs have gas-rich host galaxies. We notice that 
from all galaxies in our sample the QSO, OQ 208, is the galaxy 
with the highe st molecular ga s mass (1.4 x 10 10 M Q ). Finally 
the sample of ISolomon et al.l (1 19971) has the highest average 
of molecular gas mass of 10 10 M Q , with a sample of ULIRGs, 
where intense star formation is happening. Table [12] gives a 
summary of the average molecular gas masses compared with 
other samples, and the characteristics of the samples. 
Note that when we compare our sample with the others, we 
use the median value for the molecular gas mass that was 
calculated using only the detected galaxies for consistency with 
the other studies where the upper limits were not used for their 
calculations. 

If we compare now with a subsample of galaxies from TANGO 
itself, that have been detected in both, IRAS (for 60 ymi and 100 



Sample 


Molecular gas disk 


Dispersion 




\ 1°) 


\ lo) 


This sample 


27.6 


9.4 


Wiklind et al. (1995) 


25 


25 


Mazzarella etal. (19931 


50 


35 


Evans et al. (2005) 


55.6 


25 


Bertram etal. (2007) 


42.3 


13 



Table 13. Percentage of the galaxies, within the detected galax- 
ies of each sample of the double horn CO profile feature 



pmi) and CO(1-0); we have a total of 15 galaxies with a median 
value of 1.22xlO 9 M . This is an order of magnitude higher 
than with the median value of the molecular gas mass of the 
complete sample. It is in complete agreement with the statement 
that galaxies with higher FIR fluxes have more molecular gas 
mass than the galaxies dimmer in the FIR emission. 



6.3. Molecular gas disk 

The double-horned CO profile was previously noticed by other 
authors. We compare the fraction of double horn CO profile 
in our sample with the fraction in other samples. For instance, 
Evans et al. (2005) noted a double horn CO emission line profile 
in many of the spectra. Out of nine detected galaxies, 5 present 
the double horn profile feature. They studied two cases for which 
they had interferometric data, finding that in one case this was 
caused by a CO absorption f eature at the systemic velocity of 
the galaxy (lEvans et al.| [l999) and in the other case, 3CR 31, the 
molecular gas is distr i buted i n a molecular gas disk (Okud a et al.l 
12005 ). Wiklind et al] d 19951) ha d one case of double horn pro- 
file out of 4 galaxies detected. IBertram et al. I d2007h found 11 
double horn prof il e featu re within their 26 detected galaxies and 
IMazzarella et alJ d 19931) has 2 out of 4 galaxies detected with 
the double horn profile. In Table [13] we give the percentage of 
molecular gas disk traced by the double horn profile for each 
sample. Like for the molecular gas propert ies, the closest sam- 
ple to our radio galaxy sample is the one of Wik lind et al.l {l995) 
composed by genuine elliptical galaxies. Thus it seems that the 
presence or not of a radio source is independent of the molecu- 
lar gas properties, at the kpc scale, of the host elliptical galaxy. 
The FIR-selected samples have higher percentage of molecular 
gas disk than the sample selected only on the basis of the radio 
continuum properties. The dispersion value of the percentage of 
the molecular gas disks in each sample was computed using a 
Poissonian statistics. 

We want to add that we do not discard the possibility of the dou- 
ble horn profile in these galaxies might be caused by the pro- 
jection effect (i.e. inclination of the galaxy onto the sky plane). 



6.4. Dust 

We also compare the dust properties in our sample of radio 
galaxies, hosted by e lliptical galax i es, wi th the sample of el- 
liptical galaxies from Wiklind et al. ( 1995) (see Fig. 151 and [Til. 
The dust component i n our sample is sligh tly hotter than in the 
comparison sample of Wiklind et al. ( 1995) with a median value 
of 34 K. The media n value of the gas-to-dust mass ratio for 
Wikl ind et alJ (1 19951) is 700, much higher th an the median value 
of the ratio in our sample, which is near 260. lEvans etal] d2005l) 
found a gas-to-dust mass ratio even higher with a median value 
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Criteria of selection 


Sample Name 


Molecular gas mass 






(xlO 8 M„,) 


F'llintirnl 

JLjl 11 If \ IV. tl 1 wUldAlv.l 


This samnlp 
i ^iii 1 1 i- 1 *~ 


3.7 




Wiklind et al. (1995) 


6.8 


FIR galaxies 


Evans et al. (2005) 


37 


and QSO's 


Mazzarella et al. (1993) 


11 




Bertram et al. (2007) 


30 


ULIRGS 


Solomon et al. (1997) 


100 



Table 12. Comparison samples and the median value of their molecular gas masses 



3 -°>o- 



a) b) 

Fig. 11. Histograms of the dust mass and the dust temperature of the sample of Wikli nd et alj (Q~995). Fig. a) is the distribution of 
the logarithmic mass and Fig. b) is the distribution of the temeprature. 



of 860 . The sample of low luminosi ty QSOs (fromlBertram et al. 
120071) has a lower ratio of 452 and Ma zzarella et all d 19931) with 
their sample of radio galaxies detected by IRAS has a higher 
gas-to-dust mass ratio of 506. 



6.5. CO vs. FIR 



The sample of ISolomon et al.1 dl997l) . composed of ULIRGs, is 
the sample with the hig hest Lfir/L' cp o f 146 L Q (Kkm/spc 2 y l 
meaning, according to lSolomon & Saga (Il988l) . that this group 
of galaxies are strongly interacting/merging galaxies with a 
strong SF, although this is mainly for spiral galaxies and we 
can not use this study to classify the rest of the galaxies here. 
Our sample has a ratio of 52, which is quite diff erent from the 
ratio o f 39 in the el l iptic al galaxi es observed by IWiklind et al.1 
(119951) . lEvans etal.1 (120051) and IWiklind et al l (1 19951) seem to 
have similar ratios since Eva ns et al .1 (120051) has a ratio of 35. 
Then iBertram et al.l (120071) have a ratio of 49.1 and finally 
Mazzarella et al.l d 19931) with a very high ratio of 82. 
On Table [TT] we give the median values of the CO(1-0) and 
FIR luminosity for the comparison samples together with the 
median value of the star fo rmation efficiency (SFE -LfjR/L' C0 ). 
As alread y known, (see e g.lGao & Solomonll2004 . the ULIRG 
sample of ISolomon et al.l d 19971) has the largest SFE indicator 
with a value of 147. With a me dian value of LfiR of 1 .9 4 x 10 11 , 
the sample of radio galaxies of Mazzar ella et al.l (1 19931) exhibits 
a high SFR tracing probably the star formation in merger 
galaxies. Similarly to the ULIRGs they show a quite large 
SFE indicator (82). The sample of low luminosity QSO host 
galaxies from Bertram et al. (2007|) as well as the galaxies in 
this sample, were chosen without any bias towards the FIR 




This sample 
Evans et al. 2005 
Solomon et al. 1997 
Bertram et al. 2007 
Mazzarella et al 1993 
Young et al 95-96 
Braine et al 92 
Wiklind et al 95 
Linear Fit 
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Fig. 12. LfiR vs. L' co of all the samples used in this paper. As 
can be appreciated from this plot, galaxies from Solomon et al. 
(fl997l) which are ULIRGs, are the ones that better fit to the linear 
fit, as expected from the information in Table [TT] 



emission and their SFE appear to be similar with a value of 49 
and 52 respectively. Fina lly the sample of elliptical galaxies 
( Wiklind & Henkel 1995) and of FIR-selected radio galaxies 
( Evan s et al.ll2005l) share a low value for the SFE indicator with 
39 and 35. 
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The CO versus FIR luminosity is used in spiral galaxies to relate 
a tracer of the molecular gas, the CO emission, with a tracer of 
the star formation, the FIR emission. We compared the LfjR vs. 
L' c0 relationship in our sample with the one in other samples. 
A linear regression was fitted to the whole set of data (see 
Fig. [T2l and to the individual samples. The linear regression is 
then log(Lco)=rn-log(L/r/R)+b; where m=0.9 and b=-0.9 with a 
mean square error of 0.43. Considering the samples all together 
we obtain a;^ 2 =0.022. 

Since our data are well fitted by the CO-FIR relationship 
compared to the other samples, we can argue that the radio 
galaxies hosted by the elliptical galaxies should exhibit at least 
a low level of star formation. 

Besides the samples mentioned befo re, in Figure FPU w e in- 
clude as we l l with 3 more samples: lYoung et al.l ( 1 1995b and 
lYoung et al.l ([1996), a sample of spiral galaxies observed as 
part of the FCRAO Extragalactic CO surve y. We included 
also the sample of iBraine & Combesl dl992l) . composed of 
spiral galaxies. These samples were used to add well-defined 
samples for the CO-FIR relationship since the star formation is 
completely dominating that relationship in spiral galaxies. 



7. Discussion 

7.1. Origin of the molecular gas 

As we have shown, the molecular gas content of the host ellipti- 
cal galaxies of the powerful radio- AGN is low compared to spi- 
ral galaxies. The origin of that gas can either be internal (stellar 
mass loss, cooling gas from a galaxy halo of hot gas) or exter- 
nal (minor/major merger, accretion-fed from cosmological fila- 
ments, c ooling flows), when t he radio galaxy is at the center of 
a cluster. ISalome et aU f2008) provide an exhaustive discussion 
about the molecular gas in cooling flow which can be accreted by 
an elliptical galaxy. We estimate the stellar mass loss 77 by using 
the B magnitude of the elliptical host ( Athev et al. 2002), namely 
77 = O.OO78(M /1O 9 L , B ) M /yr. Fig. dshows the distribution 
of the stellar mass loss in our sample of radio galaxies detected 
in CO or with an upper limit of molecular gas mass lower than 
5 x 1O 9 M in order to constrain the stellar mass loss contribution. 
The median value of the stellar mass loss is 0.37 M .yr~'. This 
values implies a median time of 1 Gyr to refurbish the galaxy of 
molecular gas provided by the stellar mass loss wi thout taking 
into a ccount the consumption rate of the AGN (e.g. Davi d et alj 
120061) . As shown on Fig. [13] a substantial fraction of radio galax- 
ies have a filling time M(H-i)lr\ of molecular gas provided by the 
stellar mass loss larger than 1 Gyr (36%) reaching time-scales 
close to the age of the Universe. It all ows together with the ac- 
cretion rate of the radio-AGN (see e.g. iLagos et al. 2009) to dis- 
card the stellar mass loss as a mechanism sufficient to provide 
the molecular gas in these elliptical galaxies. 
In other words, we would expect the gas lost by the stars to be hot 
and to enrich the hot halo gas through stellar wind, in spheroidal 
galaxies; spherical geometry dilutes the gas, such that the density 
is never high enough to cool. The cooling time begins >10Gyr, 
so must of this gas remains hot and will not be seeing in H2/CO. 
This is not the case in spiral galaxies where the gravity and den- 
sity are stronger so that the gas can cool back to the HI-H2 phase. 
If the molecular gas in the radio galaxies was provided by the 
stellar mass loss we would expect, at a first order, the molecular 
gas mass to correlate with the blue luminosity, a proxy of the 
galaxy stellar mass. That hypothesis is discarded by the data on 
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Fig. 13. Histogram of the gas filling time for our sample. 



The external hypothesis for the origin of the molecular gas is 
more likely to explain a large part of the molec ular gas present in 
the elliptical galaxies hosting the radio-AGN. iLim et alj (2000) 
favored the minor merger scenario for the origin of the molecu- 
lar gas in 3CR 3 1 because of the smooth isophotes which do no 
show any signs of recent perturbations caused by a major merger. 
However, 3C R 3 1 is the brightest g alaxy in a cluster showing a 
cooling flow dCrawford et alJ [T999). 

The inspection of the optical images shows that several ra- 
dio galaxies are in dense environment where the minor/major 
mergers should be frequent enough to deposit molecular gas in 
the elliptical galaxies. The accretion of gas from cosmological 
filaments has bee n studied recently by numerical simulations 
(iDekel et a l. 2009) and should be an additional source of molec- 
ular gas in these galaxies. 

In the case of 3CR 31 there are sug gestions, bv [Martel et all 
( 1999) that the galaxy is interacting. iBaldi & CapettJ(fc 008) sug- 
gest th at 3C 264 has no young stellar population. lEmonts et al.l 
(2006) propose that B2 0648+27 f ormed from a major merger 
event that happened < 1.5 Gyr ago. lGarcia-Burillo et al ] (120071) 
state that the distortions in B2 01 16+31 may reflect that the disk 
is still settling after a merger or an event of gas accretion; or 
maybe the jet and the cone-line features might be interacting 
with t he disk producing the reported distortion. iGiroletti et al.l 
(2005) states for NGC 4278 that its central BH is active and able 
to produce jets, however, the lifetime of its components of <100 
yr at the present epoch, combined with the lack of large-scale 
emission, suggests that the jets are disrupted before they reach 
kilo-parsec scales. In conclusion, the origin of the gas in these 
galaxies can be multiple, in majority external, and a composite 
scenario can be invoked in each galaxy. 



8. Summary and Conclusion 

We present a survey of CO(1-0) and CO(2-l) lines for 52 radio 
galaxies in the Local Universe selected only on the basis of their 
radio continuum emission: 

- The detection rate is 38% (firm - 20/52) with 58% total 
(30/52), including tentative detections. 

- Indication of molecular gas disk (double-horn) is found in 
15% of the galaxies (27% of the detected galaxies) with con- 
firmation in 3CR31 with PdBI CO(1-0) map. 
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Fig. 14. Histogram of the stellar mass loss in M /yr for our sam- 
ple of radio galaxies. 

- The CO(2-l)-to-CO(l-0) line ratio is 2.3 without taking into 
account beam dilution effect, and likely to be ~ 0.6 with 
corrections, suggesting optically thick and sub-thermal CO 
emission. 

- The low ratio of /100//60 also suggest that this sample may 
be forming stars. 

We compare the properties of our sample with different samples 
of elliptical galaxies, QSO and ULIRGs: 

- We find low molecular gas content compared to the other 
selected samples, with a median value of 2.2 x 10 8 M Q . 

- We confirm that FR-II type galaxies are characterized by 
lower CO dete ction rates than FR -I and FR-c radio sources 
as claimed by lEvans et aD J2005), and that the FR-II have 
a higher molecular gas mass, as predicted from the existent 
Malmquist bias in this sample. 

- The Lfir vs. L' co relation is similar to what is found in the 
other samples suggesting that some star formation could be 
taking place in these elliptical radio galaxies. 
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Appendix A: Individual Galaxies 

We present in this appendix various properties on the 29 detected 
galaxies, with increasing 3C numbers, then NGC and B2 num- 
bers: 

- 3C 31: NGC 383, or Arp 33 1 . For information on this galaxy, 
see § 1341 

- 3C 66b: UGC 1841. The well kno wn bright r adio j et of 
this galaxy was first discovered by iNorthovei (1 19731) an d 
the optical jet was first detected by iButcher et al.l ( 1980). 
Meng & Zhoul (120061) were able to perform a detailed study 
of the jet using the HST/WFPC2 images comparing them 
with the high resolution radio images and finding excellent 
correspondence. The optical image reveals the presence of 
a disk, which is in agreement with the CO spectra of the 
galaxy. The detection is a very faint one. The central disk di- 
ameter is about 10" in the optical image, and the CO velocity 
width is 250 km/s. 

- 3C 83.1: NGC 1265 (in Abell 426). It is one of the seven 
cases where the molecular gas has been clearly detected in 
the CO(2-l) transition and not in the CO(l-O) transition. The 
o ptical image shows a dust lane of 2.24" length, according 
to iMartel etall (1 19991) . This dust lane is oriented at -171°, 
nearly orthogonal to the radio jet. If the CO emission corre- 
sponds to the dust lane, its small size explains the CO(l-O) 
signal dilution and non-detection. 

- 3C 88: UGC 2748. It is clearly detected in the CO(l-O) tran- 
sition. The optical image shows no peculiar feature. 

- 3C 129: It is another example of CO(2-l) detection with an 
upper limit in the CO(l-O) line. The FWHM velocity width 
of this galaxy is 200 km/s. There is no hint of interaction in 
the Hubble image. 

- 3C 264: NGC 3862 (in Abell 1367). This galaxy has a dou- 
ble horn profile, in both transition lines, typical of a molec- 
ular gas disk (see Figure ??). The velocity width of this 
galaxy is about 200 km/s. The disk profil e can be seen in th e 
optical image, wh ich reveals also a jet (ICrane et al ] [19931) . 



Baum et alX(l 19971) suggest that the optical synchrotron emis- 
sion, clearly visible in the optical i mage from the HST , is 
associated with the jet. According to lMartel et al.l (1 1999b the 
nucleus is unresolved, and the host galaxy projected image 
is very circular and smooth. 

- 3C 305: IC 1065. In the CO(1-0) spectra, the detection is 
very clear, with a broad velocity width of ~600 km/s. We 
have no data for the CO(2-l) e mission line. The HST image 
presented bv lMartel et al.l (1 19991) shows filamentary and dis- 
turbed swaths of dust stretching across the western side of 
the g alaxy and a twisting "arm" of emission in the eastern 
side. lJackson et al.l (1 1 995 ) noticed a 1 "scale cone-line exten- 
sion of the continuum emission north of the core, in a direc- 
tion almost 90° from the radio axis. They proposed that this 
might be an effect of obscuration of near-nuclear emission 
in all other directions rather than an example of scattering of 
anisotropic directed radiation from a hidden active nucleus. 

- 3C 321: The HST optical image shows an extended dust ab- 
sorption in the galaxy. It is clearly detected in the CO(l- 
0) transition with a velocity width of ~500 km/s. For the 
CO(2-l) emission we have no data. The HST image reveals 
a nearby galaxy that might be interacting with this galaxy. 

- 3C 327: It is clearly detected in the CO(1-0) transition line, 
with a FWHM of ~200 km/s and with one of the highest dust 
temperature from this sample (64K). 

- 3C 353: It is clearly detected in the CO( l-O) transiti o n line , 
with a velocity width of -200 km/s. IMartel et af] d 19991) 



propose that the outer isophotes of 3C 353 are very cir- 
cular while the inner isophotes are elongated in a south- 
east to northeast di rection and are roughly peanut-shaped. 
IMartel et al.l dT999) also stated that this may result from 
a small-scale dust lane bifurcating the nucleus in a rough 
north-south direction or from a true double nucleus. 

- 3C 386: It is clearly detected in CO(2-l) with a veloc- 
ity width of -175 km/s, but not in CO(1-0). The HST op- 
tical image shows a bright optical nucleus where strong 
diffraction spikes d ominate the core of the elliptical galaxy 
(IMartel et al.lll999h . 

- 3C 403: It is clearly detected in the only transition observed, 
CO(1-0), with a velocity width of ~500 km/s. It presents the 
double ho rn profile already m entioned before and shown in 
Figure ??. Martel et al. ( 1999) suggest this galaxy consists of 
two systems: a central elliptical region surrounded by a low- 
surface brightness halo with a sharp boundary at a distance 
of 3 Kpc northwest of the nucleus. In the northwest region 
of the halo, two or three very weak dust lanes are barely dis- 
cernible. 

- 3C 442: UGC 11958, Arp 169. Detected in the CO(2-l) 
emission line, but not in the CO(1-0). The HST optical im- 
age shows the elliptical shape of the galaxy but no evidence 
of dust obscuration. The isophotes of the outer halo of this 
elliptical galaxy are relatively smooth but within the central 
520 pc, they are irregular and the light distribution becomes 
non uniform. 

- 3C 449: It has a tentative double horn profile that can be 
noticed from the CO(1-0) spectrum (Figure ??), but one 
side is stronger than the other one. From the HST image of 
this galaxy it is possible to see the dust absorbi ng the visi- 
ble lig ht, not completely edge on. According to Mart el et al.l 
(1999) the morphology of this galaxy suggests that we are 
viewing the near side of an inclined, geometrically thick 
torus or disk. The velocity width is about 500 km/s. 

- NGC 315: It is detected in CO(1-0) only, there is no visible 
feature in the HST image. 

- NGC 326: Detected with a higher intensity in CO(2-l) than 
inCO(l-O). 

- NGC 541: Arp 133 (in Abell 194). This galaxy has been 
detected in CO(2-l)tmt not in CO(1-0). From the HST image 
there is another galaxy that might be falling into NGC541 
and there is another larger galaxy pr obably having some tida l 
effects on NGC 541. According to iNoel-Storr et alj d2003l) . 
this cD SO galaxy has a radio core on VLBA scales and a 
core-jet morphology on VLA scales. The central isophotes, 
measured from the WFPC/2 images, vary considerably. The 
gas does not exhibit regular rotation profile. 

- NGC 708: in Abell 262. Detected in CO(1-0) but not in 
CO(2-l). This is the central galaxy of a cooling flow clus- 
ter, detected bv lSalome & Combesl d2003l) . 

- NGC 3801: It shows a double horn profile, detected in the 
CO(1-0) spectral line, with a velocity width of ~ 600 km/s. 

- NGC 4278: Also detected bv lCombes eTall (120071) . The line 
emission is stronger in CO(2-l) than in CO(1-0), with a 
broader width. Note that t his is the closest rad io galaxy of all 
(9 Mpc). According to ICombes et"aT] {2007) the dust mor- 
phology across the disk consist in irregular patches or lanes; 
this suggests th at this galaxy has recently accreted its gas 
(ISarzi et alJl20"06l) . 

- NGC 5127: There is a clear detection in the CO(2-l) line. 
The dust absorption is not clearly visible from the HST even 
though the galaxy is not so remote (z=0.016). 
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- NGC 7052: It exhibits a double horn profile very clear in the 
CO(l-O) emission (Fig. ??), but probably also in the CO(2- 
1) emission, although it i s only a tentative detection there. 
According to iNieto et all (1 19901) the size of the dust disk is 
about 4". 

- B2 0116+31: This is a peculiar galaxy in this sample, show- 
ing a strong absorption line in CO(l-O). The molecular/dusty 
disk h as been studied in details by Garcfa-Burillo et al. 
d2007l) using the IRAM PdBI. The absorption line is sur- 
rounded by emission at both blue and red-shifted wings. 
Because of the absorption, we cannot discard a possible dou- 
ble horn profile. As mentioned in 33.21 this absorption is the 
signature of molecular gas mass on the line of sight towards 
the AGN covering a very small area and because of its small 
filling factor, this absorption is quite rare, contrary to what 
could be expected for a radio-selected sample. Its double 
horn profile is visible as well in the CO(2-l) transition line, 
where the absorption is much weaker, as well as the con- 
tinuum at this frequency, which was not detected; contrary 
to the CO(l-O) with a strong continuum flux of 164 mJy.. 
The absorption line in this galaxy should cause an underes- 
timation of the CO( 1 -0) integrated intensity and therefore an 
overestimate of the line ratio. 

- B2 0648+27: It is detected in CO( l-O) &CO(2-l) with a line 
ratio of 2.82 -the highest of all-. lEmontsi J200J&) Observed 
this galaxy with the VLA-C, where they detected HI in emis- 
sion and in absorption noticing that it is distributed in a large 
and massive ring-like structure. This distribution and kine- 
matic of the HI gas, toget her with the presenc e of faint tails in 
deep optical imaging ( H eisler & Va der 1994) and the detec- 
tion of galaxy-scale post starburst young stellar population, 
imply that this galaxy formed from a major merger event that 
happened >1.5Gyr ago. 

- B2 0836+29B: It is clearly detected in the CO(1-0) emis- 
sion line, but not in CO(2-l). Nothing peculiar in the optical 
image. 

- B2 0924+30: It shows a peculiar perturbed dust lane in the 
optical image which could possibly be due to interactions. 
It has not been detected in CO(2-l), although there is a clear 
detection of the CO(1-0) line. The HST image shows another 
smaller galaxy that might have fallen into B2 0924+30. 

- B2 1347+28: in Abell 1800. Detected in the CO(2-l) line 
but not in CO(1-0). The disk in the optical image does not 
seem to be perfectly elliptical, may be due to interactions. 

- OQ 208: Mrk 668. It has only been observed in the CO(1-0) 
line and it presents a clear double horn profile as shown in 
Figure ?? with a velocity width of about 400 km/s. 



Appendix B: Beam/source coupling 

To correct the observed emission line br i ghtnes s temperatures 
for beam dilution we know that dThi et alj (120041) 1: 



T - 

1 MB 



T MB^Qsource + ^beani) 



(B.l) 



where T MB is the real main beam temperature and Tmb is the 
observed main beam temperature. By assuming an axisymmetric 
source and beam distribution we can represent their distribution 
by an angular 6 parameter leading to a correction factor for the 
CO(2-l)-to-CO(l-0) line ratio of: 



Fig. B.l. Correction factor K for the CO(2-l)-to=-CO(l-0) line 
ration as a function of the size 6 of the source detected in CO( 1 - 
0), to take into account the beam dilution. The possible sources 
go from a point like source to a source with a size of 22" 



(see lde RiickeetalJ (120061) for more details) where 9b-230 is the 
beam size at 230 GHz and 6b-\\5 is the beam size at 115 GHz; 
9 S is the angular size of the source. We recall that for a gaus- 
sian beam D. beam = f P a a>doj = ^n9 2 h » 1.133^, although 



K = 



f- +6 2 

Vll5 ^ u s 



(B.2) 



the 1.133 cancels out in Eq. IB.2l The correction factor K for the 
CO(2-l)-to-CO(l-0) line ratio is shown on Fig. IB. II as a func- 
tion of the angular source size 6, it varies between 0.25 for a 
point source up to 1 . for a completely extended source. 
Besides for the galaxy 3CR 31, we did not apply the beam di- 
lution correction because we dont have the molecular size of 
the galaxies and therefore we do not know which correction we 
should apply. 



Appendix C: 3CR 31 

This galaxy has been observed in more details than the rest 
of the galaxies. It is an FR-I radio galaxy at a distance of 71 
Mpc hosted by an elliptical D galaxy in the Zwicky cluster 
0107.5+3212. A dou ble horn profile of the CO(1-0) and CO(2- 
1) lines was found bv lLim et al.l (2000) giving a first indication 
of a central molecular gas disk (see Fig. IC31 >. The Nobeyama 
Mi llimeter Array (N MA) CO(1-0) interferometer observations 
bv lOkuda et al.l (120051) confirmed the presence of this molecular 
gas disk. Our PdBI observations are about the same spatial reso- 
lution than the NMA ones but with a higher sensitivity as shown 
by the CO(1-0) integrated intensity map on Fig. IC.2I The mean 
surface density of the molecular gas in the center is found to be 
34OM .pc~ 2 given the spatial extension of 8"for the molecular 
gas disk as shown on the CO(1-0) intensity profile along the ma- 
jor axis (see Fig. lC.U . 

From the PdBI image (see Fig. IC.ll and IC2I > we know 
that the molecular gas extension in 3CR 31 is about 8"large. 
Applying the correction factor for the beam dilution (see Fig. 
lB5T i to the CO(2-l)-to-CO(l-0) line ratio of 2.47 found with 
the IRAM-30m observations, the actual line ratio would go 
d own to 0.8 assuming that the CO(2-l) is 8"large. According 
to iBraine & Combes! d 19921) a line ratio of about 0.7 implies an 
optically thick gas with an excitation temperature of about 7K. 
Figure IC.4I shows the position-velocity (PV) diagram for the 
CO(1-0) emission in 3CR 31. We can see that at 1" the maxi- 
mum velocity Vma.,^190 km/s. Applying the correction for the 
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Angular Offset 

Fig. C.l. CO(l-O) intensity profile along the major axis for the 
3CR 31 galaxy. 




Fig. C.l. CO( 1 -0) map of the PdBI for the 3CR 31 galaxy. 

inclination (i=39°) the maximum rotation velocity at a radius 
of l"is V rotnmx -V max / sin(/) = 317 km s _1 . Thus the dynamical 

mass inside a radius of 1 "is estimated using M f /„, = R *^ m = 
1.02 x 10 10 M Q . We estimate the total molecular gas mass in- 
side a radius of 1" to be 1.04 x 1O 8 M using the mean surface 
gas density, therefore the molecular gas represents about 1% of 
the dy namical mass in the very center of 3CR 3 1 . Okud a et alj 
(120051) discussed in detail about the dynamical and stability im- 
plications for the molecular gas in the center of 3CR 3 1 . 
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Fig.C.3. 3CR 31 spectra of the CO(1-0) (upper image) and 
CO(2-l) (lower image) lines. 
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Fig. C.4. PV-diagram along the major axis for the CO(1-0) emis- 
sion in the 3CR 3 1 galaxy from the PdBI data. 



